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Abstract 



The elliptical flow of fragments is studied for different systems at incident energies between 50 

in : 

and 1000 MeV/nucleon using Isospin-dependent Quantum Molecular Dynamics (IQMD) Model. 

OO 

Our findings reveal that elliptical flow shows a transition from positive (in-plane) to negative 
■ (out-of-plane) value in the mid-rapidity region at certain incident energy, known as transition 

o : 

energy. This transition energy is found to depend on the model ingredients, size of the fragments, 
composite mass of the reacting system as well as on the impact parameter of the reaction. A 

x 

reasonable agreement is observed for the excitation function of elliptical flow between the data 
and our calculations. Interestingly, the transition energy is found to exhibit a power law mass 
dependence. 
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I. INTRODUCTION 



The information about the nature of equation of state is still one of the burning topic 
of present day nuclear physics research in general and heavy-ion collisions in particular. A 
quite good progress has been made in the recent years in determining the nuclear equation 
of state from heavy-ion reactions l|, |2j. Among different observables, collective flow enjoys 
a special status. This is due to its sensitive response to the model ingredients that define 
equation of state. A lot of theoretical and experimental efforts have been made in study- 



ing the collective flow in heavy-ion collisions 



-|9|. This collective motion of the particles 



in heavy-ion collision can be studied via directed and elliptical flows. The directed flow, 
which measures the collective motion of the particles in the reaction plane, has been studied 
extensively at BEVALAC, SIS and AGS energies 10j. This flow is reported to diminish 



at higher incident energies due to the large beam rapidity. Therefore, elliptical flow is 
much more suited at these incident energies. The elliptical flow describes the eccentricity of 
an ellipse like distribution. Quantitatively, it is the difference between the major and minor 
axis. The orientation of the major axis is confined to azimuthal angle <j) or 0+| for ellipse 
like distribution. The major axis lies within the reaction plane for 0; while 0+f indicates 
that the orientation of the ellipse is perpendicular to the reaction plane, which is the case 
for squeeze out flow and may be expected at mid rapidity [12j. Therefore, the elliptical flow 
is defined by the second order Fourier coefficient from the azimuthal distribution of detected 
particles at mid rapidity. Mathematically, 

= Po (l + 2v 1 Cos(j) + 2v 2 Cos2(j)). (1) 

d(p 

Here <fi is the azimuthal angle between the transverse momentum of the particle and 
reaction plane. The positive value of the elliptical flow < Cos2<p > reflects an in-plane 
emission, whereas, out-of plane emission is reflected by its negative value. The reason 
for the anisotropic flow is orthogonal asymmetry in the configuration space (non-central 
collisions) and re-scattering. In the case of elliptical flow, the initial "ellipticity" of the 
overlap zone is usually characterized by a quantity e = ^ 2 ~^ 2 ^ , assuming the reaction 
plane being xz. As the system expands, spatial anisotropy decreases. From the above 
discussion, it is clear that the second order flow (elliptical flow) is better candidate for 
determining the nuclear equation of state compared to first order sideward flow (directed 
flow). 



In recent years, several experimental groups have measure the elliptical flow. The FOPI, 
INDRA and PLASTIC BALL collaborations are actively involved in measuring the 

excitation function of elliptical flow from Fermi energies to relativistic energies. In most of 



these studies, igAu 1 ^ 7 + 79AM 197 reaction has been taken 0, 5|. Interestingly, a change in 
the elliptical flow was reported from positive to negative values around 100 MeV/nucleon. 
Both the mean field and two-body binary collisions play an important role in this energy 
domain. The mean field is supposed to play a dominant role at low incident energies. The 
binary collisions starts dominating the physics gradually. A detailed study of the excitation 
function of elliptical flow in entire energy region can provide a useful information about the 
nucleon-nucleon interactions related to the nuclear equation of state. 

As discussed above, lots of attempts have already been made in the literature to explore 
different aspects of directed sideward flow. In this paper, we attempt to study the different 
aspects of elliptical flow v 2 - 

For the present study, Isospin-dependent Quantum Molecular Dynamics (IQMD) model is 
used to generate the phase space of nucleons. The article is organized as follow: we discuss 
the model briefly in section-II. The results are discussed in section-Ill and we summarize 
the results in section-IV. 



II. ISOSPIN-DEPENDENT QUANTUM MOLECULAR DYNAMICS (IQMD) 
MODEL 

n 

The isospin-dependent quantum molecular dynamics (IQMD)[13j model treats different 



charge states of nucleons, deltas and pions explicitly 13], as inherited from the VUU model 



14| . The IQMD model has been used successfully for the analysis of large number of 



observables from low to relativistic energies. The isospin degree of freedom enters into 
the calculations via symmetry potential, cross-sections and Coulomb interaction [ijj]. The 
details about the elastic and inelastic cross-sections for proton-proton and neutron-neutron 



collisions can be found in Ref. 13] 



In this model, baryons are represented by Gaussian-shaped density distributions 

fi(f,M = 4r2 ■ e- (r ^ 1(t))2 ^ ■ e-^» 2 i. (2) 
7T 2 n z 
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Nucleons are initialized in a sphere with radius R = 1.12 A 1 / 3 fm, in accordance with the 
liquid drop model. Each nucleon occupies a volume of h 3 , so that phase space is uniformly 
filled. The initial momenta are randomly chosen between and Fermi momentum(p^). The 
nucleons of target and projectile interact via two and three-body Skyrme forces, Yukawa 
potential, Coloumb interactions and momentum-dependent interactions. In addition to the 
use of explicit charge states of all baryons and mesons a symmetry potential between protons 
and neutrons corresponding to the Bethe- Weizsacker mass formula has been included. 
The hadrons propagate using Hamilton equations of motion: 

dfi d( H ) dpi d( H ) 



dt dpi ' dt dfl 

with 



(3) 



(H) = (T)+(V) 

= E^+EE / mp,t)vH?,?) 

i i j>i 

x fj (f , p 1 , t)dfdr'dpdp'. (4) 

The baryon-baryon potential V l \ in the above relation, reads as: 

yv ( r _ r) — V s J kyrme + Vy ukawa + V c J oul + Vj di + V % s 3 ym 
t x b(? - f) + 1 2 5{7 - f) P ^- 1 1 '' 



t exp(\f -f\/fi) ZjZje 2 
(If — f\/fj) |f — r| 

+t 4 \n%(p/-p) 2 + l]5(f -f) 

+tfi -jpis(n'-f-). (5) 

Qo 

Here Z { and Zj denote the charges of i th and j th baryon, and T|, T| are their respective 
T 3 components (i.e. 1/2 for protons and -1/2 for neutrons). Meson potential consists of 

Coulomb interaction only. The parameters fi and t±, , t$ are adjusted to the real part of 

the nucleonic optical potential. For the density dependence of nucleon optical potential, 
standard Skyrme-type parameterization is employed. The momentum dependence V^ di of 
the N-N interactions, which may optionally be used in IQMD, is fitted to experimental 
data in the real part of nucleon optical potential. The choice of equation of state (or 
compressibility) is still controversial one. Many studies advocate softer matter, whereas, 
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much more believe the matter to be harder in nature Il5| . As noted [16j, elliptical 
flow is unaffected by the choice of equation of state. For the present analysis, a hard (H) 
and hard momentum dependent (HMD) equation of state, has been employed alongwith 
standard energy dependent cross-section. 



III. RESULTS AND DISCUSSION 

We, here, perform a complete systematic study for the mass range between 80 and 
394 units and over full range of the impact parameter. We here simulate the reactions 
of 20 Ca 40 + 20 Ca 40 , 28 M 58 + 28 A* 58 , 41 A6 93 + 41 A6 93 , 54 Ae 131 + ^Xe 131 and 
7gAw 197 + jgAu 197 at incident energies between 50 and 1000 MeV/nucleon. In addition, the 
reactions of ^Zr m + 4 oZr 96 and 44 -Rw 96 + 44 i?w 96 are also simulated to check the isospin 



effects explicitly. As noted in Ref. 13], the relativistic effects do not play role at these inci- 



dent energies and the intensity of sub-threshold particle production is very small. The phase 
space generated by the IQMD model has been analyzed using the minimum spanning tree 
(MST) 0, [l8] method. The MST method binds two nucleons in a fragment if their distance 



is less than 4 fm. In recent years, several improvements have also been suggested [19j]. One 
of the improvements is to also imply momentum cut of the order of Fermi momentum. This 
method is dubbed as MSTM method. The entire calculations are performed at t — 200 



is 



fm/c. This time is chosen by keeping in view the saturation of the collective flow 
The elliptical flow is defined as the average difference between the square of the x and y 
components of the particle's transverse momentum. Mathematically, it can be written as: 

^=<4^4>> ( 6 ) 

Px+Py 

where p x and p y are the x and y components of the momentum. The p x is in the reaction 
plane, while, p y is perpendicular to the reaction plane. 

A positive value of the elliptical flow describes the eccentricity of an ellipse like distribu- 
tion and indicates in-plane enhancement of the particle emission i.e. a rotational behavior. 
On the other hand, a negative value of t> 2 shows the squeeze out effects perpendicular to 
the reaction plane. Obviously, zero value corresponds to an isotropic distribution in the 
transverse plane. The t> 2 is generally extracted from the mid rapidity region. The parti- 
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FIG. 1: The transverse momentum dependence of the elliptical flow, summed over entire rapidity 
distribution, at b = 0.3 for different symmetric reactions at 50 (left) and 100 MeV/nucleon (right) 
respectively. The top, middle and bottom panels are representing the free nucleons (FN), light 
charged particles (LCP's) and intermediate mass fragments (IMF's), respectively. 

cles corresponding to (Y cm /Y beam > 0.1) has been denned as projectile like (PL), whereas, 
{Xc.m./Ybeam < —0.1) constitutes the target like (TL) particles. 

In FigJU the final state elliptical flow is displayed for the free particles (upper panel), light 
charged particles (LCP's) [2 < A < 4] (middle), and intermediate mass fragments (IMF's) 
[5 < A < Atot/Q] (lower panel) as a function of transverse momentum (Pt). A Gaussian 
type behavior is observed in all cases. Note that this elliptical flow is integrated over entire 
rapidity range. This Gaussian type behavior is quite similar to the one obtained by Colona 
and Toro et al. [20]. One sees that elliptical flow is positive in the whole range of P±. The 
collective rotation is one of the main mechanism to induce the positive elliptical flow |2l| . It 
is also evident from the figure that the peaks of the Gaussian shifts toward lower values of 
P t for heavier fragments. This is due to the fact that the free and light charged particles feel 
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FIG. 2: The transverse momentum dependence of the elliptical flow, summed over entire rapidity 
distribution, for LCP's at 50 (top) and 100 MeV/nucleon (bottom), respectively. The reactions 
under study having same mass number and different atomic number. The reactions are analyzed 
with MSTM algorithm. 

the mean field directly, while heavy fragments have weaker sensitivity j^ . Furthermore, 
the peak values of v-i for the free nucleons and LCP's at 50 MeV/nucleon is 0.70, 0.411, 
0.126 and 0.27, 0.20, 0.059 for the reactions of 79 Au 197 + 79 Au 197 , 54 Xe 131 + 5A Xe 131 and 
2oCa 40 + 20C0 40 , respectively and the corresponding ratios are ~ 5.0, 3.3 and 1. The mass 
ratio of these reactions is 4.93, 3.27 and 1, whereas, N/Z ratio is 1.49, 1.42 and 1. The 
t>2 ratios are in closer agreement with the system mass ratios. The results, however, are 
different at E = 100 MeV/nucleon. Note that the peak values for the free nucleon are 0.48, 
0.34, 0.134 and for LCP's numbers, are 0.132, 0.125, 0.058. Their corresponding ratios are 
rj 2.92, 2.36 and 1, indicating a clear deviation from the mass ratio. 

To further strengthen our interpretation of the estimated t>2 ratios, we display in FigfSJ 
the reactions of aqZt 96 + 4oZr 96 and 44-Rw 96 + 44-Rw 96 under the same conditions for 
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FIG. 3: The transverse momentum dependence of the elliptical flow at E = 100 MeV/nucleon for 
the reactions displayed in Figj2j The left, middle and right panels are representing target-like, mid- 
rapidity and projectile like distributions, respectively, while, top, middle and and bottom panels 
have same meaning as that of Fig. [Tj The reactions are analyzed with MST algorithm. 

LCP's. These reactions are analyzed within MST method with momentum cut. Interest- 
ingly, the N/Z effect is more visible at E = 100 MeV/nucleon, indicating that this difference 
is not due to the mass dependence alone, but is due to isospin effect also. Our findings are 



also supported by Zhang et al. 23j, where they showed that neutron-rich system exhibits 
weaker squeeze-out flow. At low incident energy (say 50 MeV/nucleon), binary collisions 
are rare, therefore isospin in the mean field does not play role. On the other hand, around 
100 MeV/nucleon, both isospin of the mean field and binary collisions contribute, making 
isospin maximum. At higher incident energies, the role of mean field reduces. This situation 
is similar to the intermediate mass fragments, where maximum value is obtained around 100 



MeV/nucleon [24]. 



To further understand the origin of this isospin effect, the transverse momentum dependence 
of elliptical flow for target-like, mid-rapidity and projectile-like distributions is displayed in 
Figj3j From the figure, we see that isospin effect originates from the mid-rapidity region or 
in other words from the participant zone. It is also clear that the isospin effects are stronger 
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FIG. 4: The transverse momentum dependence of elliptical flow for LCP's in the mid-rapidity 
region at E = 100 MeV/nucleon. The panel is exhibiting the effect of symmetry energy on the 
40-Zr 96 +40 Zr 96 reaction. 

for LCP's compared to other fragments. This is due to the fact that heavier fragments have 



weak sensitivity towards mean field 22]. 



In FigHJ we display the transverse momentum dependence of elliptical flow for LCP's in the 
mid-rapidity region with and without symmetry energy. The effect of symmetry energy is 
clearly visible in the figure. This is in agreement with the findings of Chen et al. [^j], where 
it was concluded that light clusters production acts as a probe for symmetry energy. This is 
strengthening our agreement that elliptical flow depends on the N/Z ratio or alternatively 
the isospin dependence rather than on the size of the interacting system. 

In FigJHl we display the variation of the excitation function of elliptical flow v 2 for free, 
LCP's and IMF's over entire rapidity and mid-rapidity region. The elliptical flow is found 
to become less positive (entire rapidity) or more negative (mid-rapidity) with the increase in 
the beam energy, upto a certain energy, and then again becomes more positive or less nega- 
tive. This is due to the fact that spectators move faster after the Vi has reached a minimum 
value ]5]. This energy, at which the behavior changes is found to decrease with the size of 
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FIG. 5: The variation of the elliptical flow, summed over entire transverse momentum, with beam 
energy at b = 0.3 for different symmetric reactions over entire rapidity range and at mid-rapidity 
in the left and right panels, respectively. The top, middle and bottom panels have same meaning 
as that of FigfTJ 

the fragment. It means that the flow of heavier fragments is larger compared to LCP's/free 
nucleons at all beam energies. These type of findings are also reported by different authors 



in Ref. 16] . This is true for entire rapidity region as well as for mid-rapidity region. 
The interesting phenomena of transition from in-plane to out-of-plane is observed at mid- 
rapidity region y, 



|26| . while no transition is observed when integrated over entire rapidity re- 
gion. The energy at which this transition is observed is dubbed as transition energy (ETrans) ■ 
It means that participant zone is responsible for the transition from in-plane to out-of-plane. 
That's why free particles and LCP's, which originate from the participant zone, are show- 
ing a systematic behavior with the beam energy as well as with the composite mass of the 
system. The elliptical flow for these particles is found to become more negative with the 
increase in the composite mass of the system. Heavier is the system, more is the Coulomb 
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FIG. 6: The variation of the elliptical flow, summed over entire transverse momentum, with beam 



energy at \y\ 
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< 0.1 for jgAu +79 Au reaction. Here comparison is shown with 



experimental findings of INDRA, FOPI and PLASTIC BALL Collaborations (J-li 



repulsion and more negative is the elliptical flow. This systematics of Errans with composite 
mass of the system is discussed later. 

In Fig. El we show t>2@mid rapidity (\y\ = I Vcm I < 0.1) for Z < 2 (left panel) and for pro- 

Vbeam 

tons (right panel) as a function of the incident energy. The rapidity cut is in accordance with 
the experimental findings. The theoretical results are compared with the experimental data 
extracted by INDRA, FOPI and PLASTIC BALL collaborations 4- fj. With the increase 
in the incident energy, elliptical flow v 2 changes from positive to negative values exhibiting 
a transition from the in-plane to out-of-plane emission of nucleons. This is because of the 
fact that the mean field, which contributes to the formation of a rotating compound system, 
becomes less important and the collective expansion process based on the nucleon-nucleon 
scattering starts to be predominant. This competition between the mean field and N-N 
collisions depends strongly on the effective interactions that leads to the different transition 
energy due to different equations of state. Due to repulsive nature of the momentum de- 
pendent interactions, which leads to the suppression of binary collisions, less squeeze-out is 
observed in the presence of momentum dependent interactions (HMD) compared to static 
one (H). The maximal negative value of v-i is obtained around E = 500 MeV/nucleon with 
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FIG. 7: The impact parameter dependence of the elliptical flow, summed over entire transverse 
momentum and rapidity distribution, at incident energies 50 MeV/nucleon. The top, middle and 
bottom panels are for free particles, LCP's and IMF's, respectively. 

hard (H) and hard momentum dependent (HMD) equations of state. This out-of-plane 
emission decreases again towards the higher incident energies. This happens due to faster 
movement of the spectator matter after t>2 reaches the maximal negative value j^]. This 
trend is in agreement with experimental findings. A close agreement with data is obtained 
in the presence of hard equation of state for Z < 2 particles, while, in the presence of 
momentum dependent interactions for protons. Similar results and trends have also been 



26). 



reported by Zhang et.al. in their recent communication 

The investigation of the elliptical flow with scaled impact parameter over entire rapidity 
range is displayed in FigJTl Here the top, middle and bottom panels represent the free nu- 
cleons, LCP's and IMF's. The value of the elliptical flow Vi becomes more positive with 
the impact parameter and composite mass of the system at E = 50 MeV/nucleon, while at 
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FIG. 8: The transition energies for elliptical flow at intermediate energies as a function of the 
combined mass of the system. The upper panel is for the free nucleons, while, lower panel is for 
the LCP's. 

higher energies (not shown here), it is found to become less positive (entire rapidity) or more 
negative (mid-rapidity) with composite mass of the system. This is indicating the dominance 
of the in-plane flow at low incident energies with increasing impact parameter and composite 
mass of the system. Moreover, dominance of the out-of-plane flow at higher energies with 
small impact parameter and composite mass of the system is observed. With the increase in 
the beam energy, the expansion of the compressed zone becomes more vigorous, while, with 
an increase in the impact parameter, participant zone decreases, resulting an increases in 
the spectator region indicating dominance of azimuthal anisotropy with impact parameter. 



On the other hand, it reduces with beam energy. These observations are consistent with the 
experimental findings and with other theoretical works [lf^Q, 2^]. 

Finally, we carry out the system size dependence of the elliptical flow for free nucleons and 
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LCP's. In FigJBJ we show the transition energy Errans as a function of the composite mass 
of the system for free nucleons and LCP's. From the figure, we see that the transition energy 
decreases with the composite mass of the system as well as with the size of the fragment. 
The reason for this is that the pressure produced by the Coloumb interactions increases with 
the system size. This dependence can be fitted using a power law of the kind: 

E Trans = C(AtJ t ) (7) 

The exponent r is found to be 2 times for free particles (0.67) compared to LCP's (0.35). 
This exponent is quite smaller compared the exponent of balance energy in directed flow 



261 ] . This is due to the different origin of the balance and transition energy. The balance 
energy counter balances the mean field and N-N collisions, while transition energy is due to 
the more complex effects such as expansion of the compressed zone and shadowing of the 
cold spectator matter. 



IV. CONCLUSION 



In conclusion, we have investigated the elliptical flow of fragments for different reacting 
systems at incident energies between 50 and 1000 MeV/nucleon using isospin-dependent 
quantum molecular dynamics (IQMD) model. The elliptical flow is found to show a 
transition from in-plane to out-of-plane at a certain beam energy in mid-rapidity region, 
while no such transition is observed when integrated over entire rapidity region. This 
transition energy is found to decrease with the composite mass as well as with the size of 
the fragment. The transition energy is further parametrized in term of mass power law. In 
addition, LCP's exhibit isospin effect in the mid-rapidity region. 
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